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' Abstract 

o 

(N 

, ^ . We present the scalar-tensor gravitational theory with an exponential potential in which pauli 

metric is regarded as the physical space-time metric. We show that it is essentially equivalent to 
' coupled quintessence(CQ) model. However for baryotropic fluid being radiation there are in fact no 

, coupling between dilatonic scalar field and radiation. We present the critical points for baryotropic 

CO ' fluid and investigate the properties of critical points when the baryotropic matter is specified to 

^ - ordinary matter. It is possible for all the critical points to be attractors as long as the parameters 

^ , A and {3 satisfy certain conditions. To demonstrate the attractor behaviors of these critical points, 

CN| ' We numerically plot the phase plane for each critical point. Finally with the bound on 13 from the 

■ observation and the fact that our universe is undergoing an accelerating expansion, we conclude 

, that present accelerating expansion is not the eventual stage of universe. Moreover, we numerically 

' describe the evolution of the density parameters SI and the decelerating factor q, and computer the 

present values of some cosmological parameters, which are consistent with current observational 
data. 
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' 1 Introduction 
> . 

, The evidences for the existence of dark energy have been growing in past few years. Recently the WMAP 

^ I three result [1] has dramatically shrunk the allowed volume in the parameter space. It shows that in a spatial 

flat universe the combination of WMAP and the Supernova legacy survey(SNLS) data yields a significant 
constraint on the equation of state of the dark energy w ~ — O.OTIq qq. Though the ACDM model is still an 
excellent fit to the WMAP data, it still does not exclude other alternative model for the candidate of dark 
energy. Moreover the well-known fine-tuning and coincidence problems [2] are yet unsolved in cosmological 
constant model. This motivates a wide range of theoretic studies to explain the observation: such as the 
conventional "quintessence" scalar field[3]; the k-essence field[4]; quintom model[5]; holograph dark energy [6]; 
Born-Infeld scalar or vector field theory[7]; phantom model[8] and so on. Additionally, some authors attempt 
to modify the conventional gravitational theory instead of involving the exotic matter[9]. 

In past several years the idea that dilaton field of the scalar-tensor gravitational theory as the dark 
energy has been proposed and discussed[10]. In our previous paper[ll],we have considered a dilatonic 
dark energy model, based on Weyl-scaled induced gravitational theory. In that paper, we find that when 
the dilaton field is not gravitational clustered at small scales, the effort of dilaton can not change the 
evolutionary law of baryon density perturbation, and the density perturbation can grow from z ^ 10'^ to 
z ~ 5, which guarantees the structure formation. When dilaton energy is very small compared the matter 
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energy, potential energy of dilaton field can be neglected. In this case, the solution of cosmological scale a 
has been found[12]. In recent work[13], we consider the dilaton field with positive kinetic energy and with 
negative kinetic energy and find that the coupled term between matter and dilaton can't affect the existence 
of attractor solutions. In this paper we will study the attractor properties of the dynamical system. The 
potential we choose for investigation is the exponential form for its important role in higher-order or higher- 
dimensional gravitational theories, string theories and kahiza-klein model(see the references in Rcf[I4]). 
Though the possible cosmological roles of exponential potential have been investigated elsewhere [15], here we 
will investigate its cosmological implies in our dilatonic cosmology. With the constraint from the observation 
we conclude that the present accelerating expansion is not the eventual stage of universe. 

2 Theoretical model from scalar-tensor gravitational theory 

The action of Jordan-Brans-Dicke theory is: 

s = J d^x^m - ojr"^'^^ - A(0) - Lfi^i^m (1) 

where Lfi^id is the lagrangian density of cosmic fluid, 7 is the determinant of the Jordan metric tensor 
uj is the dimensionless coupling parameter, R is the contracted Rn^- The metric sign convention is 
(— , +, +, +). The quantity is a nontrivial potential of (j) field. When h.{(j)) ^ the Eq.(l) describes the 
induced gravity, p and p is respectively the density and pressure of cosmic fluid. However it is often useful 
to write the action in terms of the conformally related Einstein metric. We introduce the dilaton field a and 
conformal transformation as follows: 

^=^e«^ (2) 

7m'^ = '^""".9/^- (3) 
where = j^Vo, = SttG. Using Eqs.(2,3), the action (1) becomes 



-I 



d^x^[^{R{g^,) + ^g^'^'d^ad^u - W{u)) + i/i„id(^)] (4) 



where L/,„,d(Vj) = ig^-^e-^-a^^/^a^V - e-2"-i/(V,) 

The transformation Eq.(2) and Eq.(3) are well defined for some w as — | < a; < 00. The conventional 
Einstein gravity limit occurs as tr ^ for an arbitrary a; or a; ^ 00 with an arbitrary a. 

The nontrivial potential of the a field, W{a) can be a metric scale form of A((^). Otherwise, one can 
start from Eq.(4), and define W{(j) as an arbitrary nontrivial potential, g^v is the pauli metric. Cho and 
Damour et.al pointed out that the pauli metric can represent the massless spin-two graviton in scalar-tensor 
gravitational theory [16]. Cho also pointed out that in the compactification of Kaluza-Klein theory, the 
physical metric must be identified as the pauli metric because of the the wrong sign of the kinetic energy 
term of the scalar field in the Jordan frame. The dilaton field appears in string theory naturally. 

By varying the action Eq.(4), one can obtain the field equations of Weyl-scaled scalar-tensor gravitational 
theory. 

Rnu - -^g^iuR = — —{[d^ad„a - -g^udpad^a] - g^^W{a) 

+ e-«-[a^v5.v - ^g^.a^^a^v] - g^ue-^^'^vm (5) 

Aa= ^^^^ - ^e-2"<^5'"'a^V9.^i' - 2ae-2«'^y(V') (6) 
A V = -Oig^^d^^d^u + e-""^^^ (7) 
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where "A" denotes the D'Alembertian. We assume that the energy- momentum tensor T^,^ of cosmic fluid is 

(8) 

where the density of energy 

P=^V'' + e-"<^F(V) (9) 

the pressure 

P=^V''-e-«'^y(V) (10) 
In FRW metric ds^ = -dt^ + a?{t){dx'^ + dy'^ + dz"^), we can obtain following equations from Eqs. (5,6,7): 

{-/ = Y(l^' + W{a) + e-'^^p) (11) 

a + 3H&+^ = l-ae-"''{p-3p) (12) 
aa 2 

p + 3H{p + p)^^a&{p + 3p) (13) 

H = ~[&' + e-^'^{p + 3p)] (14) 

where we specify the cosmic fluid p as the baryotropic matter ph with a equation of state Pb = WbPb- 
Wb = for ordinary matter, Wb = for radiation, W{a) is exponentially dependent on a as W{(t) = Woe~^'^'^. 
The effective energy density of dilaton scalar field is pa — + W{a), the effective pressure of dilaton scalar 
field is Per = ^CT'^ — W{a) and Pa — WaPa- Here we should note that in Weyl-scaled scalar-tensor gravitational 
theory the coupling of the baryotropic matter with dilaton field a is more natural. 

We can get the solution of Eq.(13) for the density of baryotropic energy: 

Pb oc e5«(i+3«'!.)<^a-^(^+""') (15) 
For matter, Wm = 0, pm oc e^"'^a~^, and for radiation, Wr = ^, Pr oc e°"^a~'^ 

3 Critical points and the attractor solution 

In this section, we investigate the global structure of the dynamical system via a phase plane analysis. 
We define 

Then Eqs. (12-14) can be written as a plane autonomous system: 

dx \f6oL . „ , 9 „ a/SA 9 „ -j 3 , ,9 _ , 

— = ^(1 - 3wb)z^ - 3a; + + Sx^ + -(1 + Wb)xz^ (17) 

:7T7 = — 7;—xz -h -—(1 + 3wb)xz - -(1 + Wb)z 3x^z -(1 Wb)z^ (19) 
aiV 4k I A 

where N — ln{a), the constraint Eq.(ll) becomes: 

+ y2 + ^2 ^ 1 (20) 

the density parameter of a field is 

na = x^+ (21) 
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the equation of state of the dilaton field a is: 



2 2 

X — y 



= ^rh (22) 

In order to investigate the expansive behavior of scale factor a, we also represent the decelerating factor: 

9 = = ^[(1 - Wb)x'' - (1 + Wb)y'' + [wb + i)] (23) 
Using Eq.(20), we rewrite Eqs. (17-19) as follows: 

^ = -3x + Xsp-y^ + lx[2x^ + (1 + w,){l - - y^)] + ^/3(1 - 2,w,){l - x" - y") (24) 



dN yy-y2^'^ + 2^^'' +^^ + '^b)^^-'^ -ynj (25) 

where (3 = a/n. It is more convenient to investigate the global properties of the dynamical system Eqs. (24,25) 
than Eqs. (17-19). Because the case y < corresponds to contracting universes, we will only consider the 

case 2/ > in the following discussion. 

We can generally find five fixed points(critical points) where dx/dN and dy/dN both equal 0(TABLE 1) 



Points 


PI 


P2 


P3 


P4 


P5 


X 


76/3(1-3106) 


1 


-1 






6(1-Wb) 


2A-(1-3ji>6)/3 


y 











(1-^)1/2 


v'/3=(3«)6-l)^+2A/3(3u.6-l)-6(t«(,2_i) 


2A-(l-3w6)/3 



Depending on the different values ofwb, A and /3, we will study the stability of these critical points. Before 
investigating the stability of critical points, we should point out that for radiation Wb = ^, pr = e°"^a~'^, 
there are no coupling between dilaton field a and radiation in Eq.(ll). Since = 3pcr, the right of Eq.(12) 
equals zero and therefore the radiation density does not contribute to Eq.(12). In this case, the form 
of Eqs. (11, 12, 24, 25) are very similar with conventional scalar field and the corresponding stability of the 
critical points are presented in[14]. Therefore we only pay attention to the case Wb = O(matter). In this case, 
we will find that the properties of critical points are more distinctive than the case Wb = ^. We present the 
critical points and its properties in following TABLE 2: 



Points 


X 


y 


Stability 


f2<x 




Q 


PI 







X- <P< V6{X > Vg) 
or -VG < (3 < X+{X < -VG) 


6 


1 


4 


P2 


1 





X > Ve and P > VG 


1 


1 


2 


P3 


-1 





X < -VG and /3 < -^6 


1 


1 


2 


P4 




(1 - f )V2 


< 6 and P'^ < 24 


1 


3 


^-1 

2 ^ 


P5 




(1 - ^)V2 


max{p-,—\/6) < 

X < min{(3+, VQ)) 


1 


^-1 
3 ^ 


1 

2 ^ 


P6 


2A-/3 


r/3^-2A/3+6il/2 
1 (2A-/3)2 i 


see Eq.(31) 


A 


B 


C 



where A = ^^ix^g^, B = -I + ^.J^^^^^ and C = Q^^^^. We define A± = A ± ^/X?^, 

P± = • 

Prom Table 2, we know that for every critical point, they are all possible to be stable, but for every value 
of the parameters A and (3, there is one and only one stable critical point; the others arc unstable. 

PI is a stable point (Figl), which corresponds to a late-time decelerating rolling attractor solution where 
neither the dilaton nor the baryotropic fluid entirely dominates the evolution. Note that 0^ = ^ ' which is 
irrelevant to the potential parameter A. 
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The properties of P2(Fig2) and P3(Fig3) are very similar, they both correspond to a late-time decelerating 
attractor solutions. In this case the right of Eq.(ll) is dominated by the kinetic energy of the dilaton field, 
which are different with the critical point PI. All the three points describe the dilaton field behaving as a 
"stiff' matter with a equation of state = 1. These three points are obtained from j/ = in Eq.(25). 
However if _ 

'^Xx+^[2x^ + {l+Wb){l-x^-y^)]=0 (26) 
we can obtain other critical points. From Eqs.{24:,25){dx/dN = 0, dy/dN = 0) we get following equation: 

{(3 - 2\)x^ - ^ [A(/3 - 2A) - 6]a; - A = (27) 

Solving Eq.(27), we can obtain critical points P4-P6. there are two different cases: 

Case 1 If /3 = 2 A, there exists only one solution(P4, Fig4) for Eq.(27). It is a stable node point for 
|A| < VG and corresponds to an accelerating expansive universe for |A| < \/2. 

Case 2 If /3 ^ 2A, there are two solutions for Eq.(27) and therefore exist two critical points(P5, P6). 

For critical point P5 it is stable for maa;(/3_, — \/6) < A < mm(/3+, \/6). P5 corresponds to a late-time 
attractor solution(Fig5) where the density of dilaton field will dominate the universe. The universe will 
undergo a stage of accelerating expansion if |A| < \/2- 

There are two constraints on the parameters (3 and A for point P6. Since t/5 = we have 

the first constraint: 

/J^ - 2A/3 + 6 > (28) 
The second constraint is from that: 0(, = 1 — fi^ > 0. So we have: 

2A^ - A/3 - 6 > (29) 

To be a stable point for P6(Fig6), the parameters A and (3 have to satisfy another constraint (from the 
analysis of the stability of Eqs.(24,25)): 

(2A-/3)[A-i(/32-/?3)i/3_ >0 (30) 

From Eqs. (28-30) we have following constraints on A and (3: 

A>/3+ for /3<0; 

A</3- for /3>0; 

/?->A>| + | for 0> I3> -V6; 

| + f>A>/3+ for V6>f3>0; 

Eq.(31) is the condition for P5 being a stable point. 

If the parameters A and /3 also satisfy another relation: 

(2A-/3)(A + /3) <0 (32) 

the stable point will correspond to a late-time accelerating expansive universe where the dark energy can not 
dominate the universe entirely. From Eqs. (31, 32), we get the range of A and /3 for an accelerating expansive 
universe: 

-/3>A>/3+ for P<-V2 

f3->X>-(3 for f3>V2 ^ ' 

However if 2A^ — A/3 — 6 = 0( A = in Eq.(27)), the two different critical points(P5 and P6) become one 
same point. In this case we find that the largest eigenvalue for linear perturbations vanishes and we must 
consider higher-order perturbations about the critical point to determine its stability. 
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In order to clearly understand the attractor properties of these six critical points, we plot their phase 
planes when the values of A and j3 lie in the corresponding range presented in Table 2. 




Figl. The phase plane for A = 3, /3 = 2(P1). 
The late-time attractor is rolling solution with 
x\ = \/6/3, 2/1 = 0,where Oo- = 2/3, q = 3/2 and 

Wtr = 1 



property of the syslerri 




Fig2. The phase plane for A = 3, /3 = 4(P2). The 
late-time attractor is a dilaton field dominated 
solution with 0:2 = 1, t/2 = 0, where = 1, g = 2 
and Wct = 1 




properly of Ihe system 



I i J J / / / . 

I i J J / / y , 

I i J / / / 

I I J ^ / yy. 

I J / / 




■y y / / /M\\\\\ 
- ' y /■ / r t \ \ WW 



Fig3. The phase plane for A = -5, /3 = -6(P3). 
The late-time attractor is a dilaton field dominated 
solution with 0:3 = — 1, j/3 = 0, where fi^ = 1, 9 = 2 
and Wct = 1 



Fig4.The phase plane for /3 = 2A = 2(P4). The 
late-time attractor is a dilaton field dominated 
solution with xa = -^/l/G, j/4 = where 0^ = 1, 

q = —1/2 and = —2/3 
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FigS.The phase plane for A = 1/3, j3 = 8(P5). 
The late-time attractor is a dilaton field dominated 
solution with xr, — ■\/l/54, j/5 — 1^53/54, where 
= 1, g « -0.944 and Wa = -0.963 



Fig6.The phase plane for A = 3, /3 = 1(P6, deceler- 
ating expansion). The late-time attractor is scaling 
solution with X(, — a/6/5, j/e = 1/5, where Q.^ = 7/25) 
q = 4/5 and w„ = 5/7 



4 Constraint from observation 

Firstly we consider the observational data that the present solar-system gravitational experiments set 
(at the 1(7 confidence level) a tight upper bound on 0loiar-system < 10~^[17]. So we have: 

< 10"^ (34) 

Therefore with the constraint on (3 and the attractor conditions presented in Table 2, we know that there 
only exists four possible late-time attractors(Pl, P4, P5, P6), corresponding to four different cosmological 
destinies: 

A. the critical point PI, it is a stable node for A_ < /3 < \/6(A > \/6) or — \/6 < /3 < A+(A < — \/6). 
Considering the upper bound Eq.(34) on /?, it requires |A| > 94.88. In this case, the universe will be nearly 
entirely dominated by ordinary mattcr(f7c = ^ < g(^) and end with a decelerating expansion. Therefore 
the present accelerating expansion is only a transient regime. 

B. the critical point P4. in this case /? = 2A, so from Eq.(34),we have A^ < | x 10~^. In this case 
the universe will be entirely dominated by dilatonic scalar field (f^o- = l)and it is undergoing an accelerating 
expansion with the value of w — 1 and q —1. Since the current density parameter of dark energy 
VLde ~ 0.7, the current cosmological stage is not a final evolutive stage yet. 

C. the critical point P5. It is stable for max (/3_,-V6) < A < mm(/3+,V6)). From Eq.(34), the 
constraint for A is 1.73 > A > —1.73. In this case the universe will be entirely dominated by dilatonic scalar 
field(OCT = !)• Whether the universe finally undergoes an accelerating expansion is determined by the value 
of \{q — |(A^ — 2)). Clearly the current cosmological stage is not a final evolutive stage yet. 

D. the critical point P6. From Eq.(33) we know that there is no accelerating expansion if (3 satisfies 
Eq.(34), so the current stage is still not a final evolution in this case. 

The second constraint is from the fact that our universe is now undergoing an accelerating expansion and 
is dominated by the dark energy. The contribution of dark energy is ~ 0.7 and its state parameter(a;o-) 
is very close to —1. So, we can exclude the points PI and P6, which do not admit accelerating expansion. 
Therefore there only exist two critical points P4 and P5 which satisfy above constraints. So we can conclude 
that the universe will finally be completely dominated by dilatonic scalar field, hence the current cosmological 
state is not our ultimate destiny and we are all on the way toward final late-time attractor. Furthermore, 
We plot the evolution of density parameter and Or(Fig7, where we consider the baryotropic matter 
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Pb as matter and radiation, not only one). Fig8 is the evolution of the decelerating factor q with respect to 
A'^, which quantitatively indicates that the universe is speeding up at present but underwent a decelerated 
regime at the recent past. Since T^q ~ 5.64(f^o^^)eV' =i 2.843 x 10'*if(where we set h 0.71), Tq ~ 2.7K, 
o-i = o-eq = 1, the scalc factor at the present epoch oq would be nearly 1.053 x lO"', then we know at present 
epoch A^o — Inao ~ 9.262. We calculate the present value of il^,, ~ 0.701, decelerating factor qo ~ —1.085 
and the transition rcdshift Zt — 0.670 in our model. It is consistent with current observation data. 
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FigT.The evolution of density parameter Q with 
respect to A^. A = 0.01, (5 = -0.02. Solid line is for 
the dilaton field fla, dotted fine is for radiation Qr 
and dashed line is for matter flm- 
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FigS.The evolution of decelerating factor q with 
respect to A'^. It clearly shows that the universe 
underwent a decelerated regime at the recent past 
with the transition redshift Zt — 0.670 



5 Summary 

In this paper, we studied the global properties of dilaton universe in the presence of dilatonic scalar field a 
couphng to a baryotropic perfect fluid. We find the general solution( Eq.(15)) for the density of baryotropic 
energy. We show that if the baryotropic matter is radiation(?i;f, = 1/3) the stability of critical points are 
the same as the scalar field[14] which is no coupling to baryotropic matter. In this case there are in fact no 
coupling between a filed and radiation. However, if the baryotropic fluid is matter(z«;, — 0), the situation is 
quite different. Since there exists coupling between a field and matter, which is mathematically equivalent 
to the coupled quintessence model. It is possible for all the critical points to be stable but for every value 
of the parameters A and /3, there is one and only one critical point to be stable; the others are unstable. 
We plot the phase plane for each critical point and the results clearly show the attractor properties of these 
critical points. With the constraints from the observation data, we conclude that our present universe does 
not reach the ultimate cosmological state yet and we are all on the way toward final late-time attractor. 
Finally we present a numerical computations on the density parameter n and the decelerating factor q. All 
the results are consistent with current observation data. 
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